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IMPROVING THE ACCURACY OF ANGULAR -MOMENTUM PROJECTION 

by William F. Ford 
Lewis Research Center 

SUMMARY 

A connection is established between Ullah's new method of angular-momentum pro- 
jection and the conventional Hill-Wheeler method. They are studied for the case where 
series truncation of some sort is required. It is shown that for a particular choice of 
angles the analysis simplifies greatly and at the same time leads to reduced truncation 
error. 


INTRODUCTION 

Recently Ullah proposed a new method (ref. 1) for performing angular-momentum 
projection, designed to avoid the difficulties associated with numerical integration of 
rapidly oscillating rotation matrices. Considering the case of axially symmetric intrin- 
sic wave functions he used LSwdin's representation (ref. 2) of the projection oper- 
ator P to write 


-l) n C J-K+n ■ . • - (i) 

n!(2J + n + 1) ! 

where fi is some rotationally invariant operator (usually H or 1) and 

C n = <®kI £jj -J"I* K > ( 2 ) 

(Both and C_ depend on K, but for simplicity this dependence will be left implicit.) 
j n 

Rather than evaluate equation (2) directly, Ullah showed that the coefficients C n could 
be found by means of a generating function 




n j s < $ kI 0 P jkI*k^ 

t 

= (2j + 1) (j + k) ! y i 

( J - K) ! 


J -J 
max 


n=0 




and then devised an elegant transformation to reduce the left side to 


(1 - x 2 )- K 





(4) 


where 

si,, G •') • ' 

If equation (3) is considered for a set of N + 1 different values of X, a simple matrix 
inversion (ref. 3) will yield the coefficients C . 


CONNECT I ON W ITH HARMONIC ANALYS I S METHOD 


Basically, Ullah’s method is an improvement of one developed previously by 
Mihailovic, Kuhawski, and Lesjak (ref. 4). Working from an expression which can be 
written (for rotationally invariant operators) in our notation as 



6J 6 J 

e " fie 
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they obtain the coefficients by means of harmonic analysis and then numerically 
solve the set of linear equations 




J=K+n 


(J + K + n)f(J - K)f 
(J - K - n) ! (J + K) ! 




( 6 ) 


to find the desired quantities fij. This is unnecessarily cumbersome, however, since 
by writing > ■ 


2 



max 


pi _ \ (J + L + n) ! 
L-K+n ” Z.J (J - L - n)! 


(J - K)? 
n) ! (J + K)! 




J— L+n 


( 7 ) 


and using the identity 


J-L 


(2L + 1) Y ( J -t L -t. n )-! illl 

Z — t (J - L - n).! n!(2L ,+ n + 1)! 
n=0 
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( 8 ) 


we may easily solve the linear equations analytically, obtaining 


J max" L 


= (2L + 1) ( L + K ) • V W 

u (L-K)! n!(2L + n 

n=0 


+ 1 )! 


'L-K+n 


(9) 


the form proposed by Ullah. The identity (8) may be established by noting that the left 
side is proportional to the hypergeometric function 2 F i(L + J + 1, L - J; 2L + 2; 1). 
This in turn may be written in closed form, and the left side reduces to 
(2L + 1)/(L + J + 1)(L - J)!(J - L) ! 

Another difficulty is that the transformation used by Mihailovic, Kujawski, and 
Lesjak to render their equations harmonic requires evaluation of equation (5) at 2N + 1 
different complex values of 6 to obtain the N + 1 coefficients C^. Since calculation 
of the overlap integral in equation (5) is difficult and lengthy on even the fastest compu- 
ters, the direct method proposed by Ullah, which requires only N + 1 different evalua- 
tions, is to be preferred. 

There still remains a problem, however. If the dimensionality of the model space 
used to define is allowed to increase, or an "inert" core is allowed to become 
"active", the value of J max can become enormous. In that case the expansion (3) 
must be truncated for practical reasons, and it is not clear how to choose the set of X 
values so as to minimize the truncation error. 

CONNECTION WITH HILL-WHEELER METHOD 


It is instructive to compare Ullah’ s method with the conventional one employing the 
Hill-Wheeler integral (ref. 5), 
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d KK^.^ sin £ d ^ 


( 10 ) 


where 


4 k<» * (•« 


-i0J, 


JK, 


is an element of the reduced rotation matrix. It has been pointed out by Ripka (ref. 6) 
that direct numerical integration of equation (10) is likely to be inaccurate for large 
values of J, because of the rapid oscillations of d^^O). The factor 
(<f> K |fi exp(-i/3jy) |$^) , on the other hand, is relatively smooth except for a strong peak 
at 0 = 0 (and possibly at 0 = n) (refs. 7 and 8). This suggests an expansion in powers 
of some suitable function of 0, and from Ullah's work (see eqs. (3) and (4)) it is appar- 
ent that one such expansion, with only a finite number of terms, is 
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-i0J v 

fie y 
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Using equation (11) in (10) results in 


°J - 


2J + 1 
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JK 
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( 12 ) 


where 


q JK (;.l) 
(n 


^ f7sin±0 

'■ ? Jo ' 2 , ' ; 


2n 


2K< 

cos - 0 ) d^ K (0)sin 0 d0 



1 - x\" n + x\ 2K p» k 2K) (x)dx 


(13) 


Here the reduced rotation matrix has been replaced by its representation in terms of a 
Jacobi polynomial (ref. 9). Because of the orthogonality property of Jacobi polynomials, 
the integral in equation (13) vanishes unless n >: J - K. Its nonvanishing values are 
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given by (ref. 10) 
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JK 

J-K+n 


2 ( _i) n (J + K)l 1 

(J - K)! n!(2J +' n + 1)! 


and we recover Ullah's result. 


(14) 


A SIMPLIFIED METHOD BASED ON ORTHOGONAL POLYNOMIALS 


The previous section illustrates the connection between JJllah's expansion and the 
conventional Hill-Wheeler integral treatment. The orthogonality i feature, however, sug- 
gests a different expansion of the integrand factor in equation (10), namely, that. the 
terms in equation (11) be rearranged so that 
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1 \ 2K 
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- £ a j 4 k<» 

J=K 


(15) 


We then obtain the remarkably simple result 


" < $ Kl nP JK^K> " a J 


(16) 


Of course, this is really not so surprising, since equations (10) and (15) together form a 
statement of the expansion theorem for orthogonal functions, that is, one implies the 
other. The importance of equation (15), however, is not only that its coefficients 
aj = fij may be obtained by exactly the same procedure suggested by Ullah, namely, 
matrix inversion, but in addition that for a particular choice of angles it is possible to 
arrange matters so that the matrix inversion can be done trivially and yet furnish a 
superior approximation to the integral if the expansion is truncated at J = J Q . 

This is accomplished by using the roots of the first neglected rotation matrix 

Jo+1 

d KK (/ V = 0 n = 0, 1, N o = J o -K. (1.7) 

(excluding the 2K-multiple root (3 = n), because for these angles the Christoffel-Darboux 
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formula (ref. 11) for Jacobi polynomials yields an "orthogonality” relation: 


V-/ 

X i^ir) d KK< = ° 


m =/ n 


( 18 ) 


J=K 


With this relation the N Q + 1 linear equations obtained by setting /3 = 0 n in equation (15) 
are easily solved, yielding 
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Z w - 
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where 


W n = 


1-1 


o 

I 

J=K 


2 J + 1\ ,J y 2 
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( 20 ) 


Since the Christoffel-Darboux formula is of the form 


N 


z 


F (x, y) 

C n P n (x)P n (y) = -j— 


n 


it follows that if y = x + e with e — 0, then 

‘ N , • 


Z c - p " 


9 d F T.T 

<x) 2 = N 


. .. 3y 


y-x 


This may be used to sum the series in equation (20), with the result 

" •* 12 


w „t2 


(J Q + l)sin 0 n 


(J„ * 1 + K)(J„ + 1 - K)d^(g n ) 


(21) 


These formulas have the appearance of an algorithm for the numerical integration 
of equation (10) and in fact are identical to those which would be obtained from; a Gauss 
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quadrature method of order N Q based on the orthogonal functions d^ K (/3). As such, 
they are known to be exact if the integrand can be expressed in terms of the first. 2 N q + 1 
such functions. We conclude that, if the expansion (15) is terminated at J = J Q and its 
coefficients are determined by equation (19), truncation error will affect only those 
for which J > 2 J q + 1 - J max ; those with J < 2J Q + 1 - J max will remain exact. The 
magnitude of the truncation error is difficult to estimate, but the Gauss method is known 
to be remarkably accurate. 

CONCLUDING REMARK 

One caveat should be sounded: if the dimensionality or number of nucleons is large, 
the number of terms which can be taken is severely limited by the computing time re- 
quired to evaluate exp(-i/3jy) . At the same time, the behavior of the inte- 

grand places a greater burden of accuracy on the series expansion. In a typical calcula- 
20 

tion for Ne, for instance, with all nucleons active and states up to lgg /2 included, the 
left side of equation (15) falls from 1. 00 at j3 = 0° to 0. 070 at /3 = 45° to 0. 001 at 
j3 = 65°. Here the strong forward peaking is of more concern than oscillations of 
d KK^’ 311(1 11 may be P rofita ble to consider a transformation on the variable of integra- 
tion before attempting numerical procedures. In that event it would probable be best to 
treat the transformed integrand by means of a Gauss quadrature based on Chevyshev 
polynomials, because of the ease in calculating weights and abscissas. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 2, 1973, 

503-10. 
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